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Abstract 

A /3-skeleton is a proximity undirected graph whose connectivity is 
determined by the parameter ft . We study /^-skeleton automata where 
every node is a finite state machine taking two states, and updating its 
states depending on the states of adjacent automata-nodes. We allow 
automata-nodes to remember their previous states. In computational ex- 
periments we study how memory affects the global space-time dynamics 
on /3-skeleton automata. 
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1 Introduction 

In computational geometry and geometric graph theory, a /3-skeleton or beta 
skeleton is an undirected proximity graph defined from a set of points in the 
Euclidean plane. Two points p and q are connected by an edge whenever their 
(3 neighborhood is empty [24 . In lune-based /3-skeletons, the /3-ncighborhood is 
defined as : the intersection of two circles of radius d(p, q)/2(3 that pass through 
p and q, if (3 € [0, 1] ; the intersection of two circles of radius (3d(p, q)/2 centered 
at the points (I - /3/2)p + ((3/2)q and {(3/2)p + (1 - /3/2)q, if /? > 1 . 

Figure [l] shows three /3-skeletons with increasing /3 parameter value, based 
on the same ten nodes. The figure shows also the /3-neighborhoods of the node 
labeled I (upper-right). In the transition from (3—0.9 to (3—1.0, the 1-node loses 
its 1-6 and 1-8 links, whereas from (3—1.0 to (3=1.5 it loses the 1-5 link. 

/3-Skeletons belongs to a family of proximity graphs, which are monotonously 
parameterised by parameter the (3. The structure of proximity graphs repre- 
sents a wide range of natural systems and is applied in many fields of science. 
Few examples include geographical variational analysis [3TJ [37J G£I| > evolutionary 
biology , simulation of epidemics [3B] , study of percolation [TH] and magnetic 
field [35], design of ad hoc wireless networks H5J [3H [3T] HU] . Thus developing 
and analysing computational models of spatially-extended systems on proxim- 
ity graphs will shed a light onto basic mechanisms of activity propagation on 
natural systems. 

Automata on /3-skeletons were originally introduced in [5] and studied in a 
context of excitation dynamics. In present paper we develop ideas of [5] along 



Alonso-Sanz R. and Adamatzky A. On beta-skeleton automata with memory. 
Journal of Computational Science 2 (2011) 1, 57-66. 



13=0.90 



0=1.00 




K=3.i 



K =2.60 



7f=2.00 



Figure 1: Three /3-skeletons on the same ten nodes. 



lines of memory-enriched automata and global dynamics. The paper is struc- 
tured as follows. In Sect. [2] we define automata on /3-skeletons. Global dynamics 
on automata for j3 > I (planar graphs) are studied in Sects. pOj and 1 2.3] , whereas 
the /3 < 1 (non-planar graphs) is studied in Sect. 2.2. Section [4] demonstrates 
the effects of weighted memory on the behaviour of automata networks. Au- 
tomata networks with non-parity rules are briefly tackled in Sect. [5j Possible 
applications of our computational findings are discussed in Sect. [6j 



2 Automata on beta-skeletons 



In the automata on beta-skeletons studied here, each node is characterized by 
an internal state whose value belongs to a finite set. The updating of these 
states is made simultaneously (d la cellular automata) according to a common 
local transition rule involving only the neighborhood of each node [5]. Thus, if 

(T) 

a\ ' is taken to denote the state value of node i at time step T, the site values 
evolve by iteration of the mapping : c t - T+1 ' = ^{^W^} € ■Mj j where A/i is 
the set of nodes in the neighborhood of i and <f> is an arbitrary function which 
specifies the automaton rule. This article deals with two possible state values at 
each site: a € {0, 1}, and the parity rule: <j\ — y, a ) mod 2 . Despite 

its formal simplicity, the parity rule may exhibit complex behaviour [23] . 

In the Markovian approach just outlined (referred as ahistoric), the tran- 
sition function depends on the neighborhood configuration of the nodes only 
at the preceding time step. Explicit historic memory can be embedded in the 
dynamics by featuring every node by a mapping of its states in the previous 
time steps. Thus, what is here proposed is to maintain the transition function <f> 
unaltered, but make it act on the nodes featured by a trait state obtained as a 
function of their previous states: a\ ^({ s j } s - being a state 

function of the series of states of the node j up to time-step T. We will consider 
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here the most frequent state (or majority) memory implementation. Thus, with 
unlimited trailing memory : s\ T ^ — mode , of 1 \ . . . , whereas with 

(T) / (T) (T— 1) (T)\ 

memory of the last r state values : s\ — mode[<j\ ' ,a\ , . . . ,a\ ') , with 
T = max(l, T — T + l). In the case of equality in the number of time-steps that 
a node was and 1, the last state is kept, in which case memory does not really 
actuate. This lack of effect of memory explains the lower effectiveness of even 
size r-memories found in the results presented below. 

A FORTRAN code working in double precision has been implemented to 
perform computations. The core of the code is given in Table [TJ The wiring of 
the n nodes is encapsulated inanxn adjacency matrix named ADJ, so that the 
iadd variable adds the states of the nodes linked to the generic node i. Thus, 
the parity rule is implemented by means of the mod 2 operation on the iadd 
variable, generating NEW(i), i.e., a\ . Previously, the MEMORY subrutine im- 
plements the itau-majority memory computing. MEMORY provides the MODE(i) 

(T) 

trait states, i.e., , on which the iadd is calculated. 

DO it=l,maxit 

call MEMORY (MODE , NEW , ISIG , n , it , maxit , it au) 
DO i=l,n;iadd=0 

do j=l,n;if (ADJ(i,j)==0)cycle 
iadd=iadd+MODE(j) 

enddo 

NEW(i)=mod(iadd,2) 
ENDDO 
ENDDO 

subr out ine MEMORY (MODE , NEW , ISIG, n, it , max it , it au) 
integer MODE(n) ,NEW(n) ,ISIG(n, maxit) 
ISIG ( : , it ) =NEW ; M0DE=NEW 
itin=max(l , it-itau+1) ; itaux=min(it , itau) 
do 1=1, n; inn=sum(ISIG(i , itin: it) ) 

if (2*inn>itaux)M0DE(i)=l 

if (2*inn<itaux)M0DE(i)=0 

enddo 

END 

Table 1: The core of the FORTRAN progam. 
2.1 The 0=1 case. 

Figure [2] shows the initial evolving patterns of a simulation of the parity rule on 
a /3=1 skeleton (or Gabriel maps) with N = 10 2 nodes distributed at random 
in a unit square. Red squares denote node state values equal one, black squares 
denote zero state values. The effect of endowing nodes with memory of the 
majority of the last three states is shown at T—A . The encircled node exem- 
plifies the initial effect of memory. Dynamics of perturbation spreading looks 
unpredictable and quasi-chaotic due to existence of long-range connections. 

Figure [3] shows the evolution of the changing rate (the Hamming distance 
between two consecutive patterns) in eleven different (3=1 simulations based in 
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Figure 2: A simulation up to T=4 of the parity rule on a (3=1 skeleton, N = 100. 
Red squares denote node state value equal one, black squares denote zero state 
values. 



one thousand nodes distributed at random in a unit square. The red curves 
correspond to the ahistoric simulations, in which case the parity rule exhibits a 
very high level of changing rate, oscillating around 0.5 . Figure [3] shows also the 
effect on the changing rate of endowing nodes with memory of the last r state 
values (blue lines). The inertial effect of memory tends to reduce the changing 
rate compared to the ahistoric model, particularly when (3 is odd. With high 
memory charges, such as r=19 in the lower left panel, the changing rate tends to 
vary in the long term in the [0.1, 0.2] interval, after an initial almost- oscillatory 
behaviour which ceases by T > r+l=20. With unlimited trailing memory 
(lower right panel) , this oscillatory pattern is never truncated, so that a rather 
unexpected quasi-oscillatory behaviour turns out with full memory. With no 
exception, the proportion of node states having one given state value {density), 
oscillates near to 0.5 regardless of the model considered. 

Figure [4] shows the evolution of the damage rate, i.e., the relative Hamming 
distance between patterns resulting from reversing the initial state value of a 
single node, referred to as damage (or perturbation) spreading. The damage 
propagates very rapidly without memory {butterfly effect), so that by T=30 the 
red curves already oscillate around the fifty per cent of cells. When memory is 
introduced to the system, the spread of damage is depleted, albeit the restraining 
effect of memory is very low if the charge of memory is low (t=3,4). With higher 
memory lengths, e.g., r=9, the depletion in the advance of the damage becomes 
apparent, though by T=150 the damage rate also reaches the 0.5 level in every 
simulation in Fig. [4]. Similar evolution of damage is found with higher limited 
trailing memories, such as t=11 and r=13 (not shown in the figure). On the 
contrary case, unlimited trailing memory appears very effective in the control 
of damage, as shown in the lower-right panel. 

Figure [4] also shows the relative Hamming distance of the ahistoric patterns 
to the corresponding historic ones. After a very short transition period, this dis- 
tance reaches a fairly permanent plateau level, oscillating around 0.5 regardless 
of the charge of memory endowed. 



4 



Alonso-Sanz R. and Adamatzky A. On beta-skeleton automata with memory. 
Journal of Computational Science 2 (2011) 1, 57-66. 




10 20 30 40 50 60 70 80 90 100 



Figure 3: Changing rate in eleven simulations up to T=100 of the parity rule 
on j3=l skeletons, N = 10 3 . Red plots correspond to the ahistoric simulations, 
blue plots correspond to r-majority memory. 
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Figure 4: Damage spreading in the scenario of Fig. [3]. Color code: black, dis- 
tance to the ahistoric pattern; red, damage spreading with no memory; blue, 
damage spreading with memory. 



2.2 A p <1 case. 

Figure[5] shows a simulation up to T=4 of the parity rule on a /3=0.9 skeleton 
based in the same nodes and initial states as in Fig. [2]. In correspondence with 
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the lower (3 value, there are more links connecting nodes in the non-planar graph 
of Fig. [5] compared to those in Fig. [5]. 

NO MEMORY MAJORITY MEMORY 



\ ■ ■ \ ■ ■ s ■ ■ ■ s ■ ■ ■"■ ■ s " . ■". ■ 




Figure 5: A simulation up to T—4 of the parity rule on a /3=0.9 skeleton. The 
location of nodes and their initial states are those of Fig. [2] 

Figure [6] shows the evolution of the changing rate in eleven different /3=0.9 
simulations, based in the same nodes and initial states as in Fig. [3]. In this 
scenario, with high connectivity, low memory charges, e.g., r=3,4, have not 
an apparent effect on the changing rate, so that there are not presented in the 
figure. Memory charges of r=5 and t—7 have a limited effect, whereas r=9, 19 
already have an apparent effect. With unlimited trailing memory, the oscillatory 
patterns get a notable amplitude, albeit below 0.5 . 
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Figure 6: Changing rate in eleven simulations up to T=100 of the parity rule 
on /3=0.9 skeletons. N=1000 . Color code as in Fig. [3]. 

Figure [7] shows the damage spreading in the scenario of Fig. [6]. In this 
highly connected scenario, memory turns out rather ineffective in the control 
of damage. Even with unlimited trailing memory, the damage grows up to 0.5 
fairly soon in most simulations, or is just delayed up to approximately T=100 
in some others. 
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Figure 7: Damage spreading in the scenario of Fig. [6]. Color code as in Fig. [4]. 
2.3 A p >1 case. 

Figure [5] shows a simulation up to T=4 of the parity rule on a lune-based /3=2 
skeleton (also termed relative neighborhood graphs) based in the same nodes 
and initial states as in Fig. [2]. In correspondence with the higher /3 value, there 
are less links connecting nodes in Fig. [8] compared to those in Fig. [2]. 
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Figure 8: A simulation up to T=4 of the parity rule on a /3=2 skeleton, N = 10 2 . 
The location of nodes and their initial states are those of Fig. [5]. 

Figure [9] shows the evolution of the changing rate in eleven different /3=2 
simulations based in the same nodes and initial states as in Fig.[3]. In this 
scenario, even low memory charges, e.g., r=3,4, have an apparent effect on 
the changing rate, and higher ones, e.g., r=9, 19 , led this parameter to very 
low values. Thus for example with r=19 the changing rate varies in the [0,0.1] 
interval. With unlimited trailing memory (lower right panel) , the non-truncated 
oscillatory patterns soon reach a very high level near 0.5 . 

Figure [lO] shows the latest snapshots of the changing nodes in one of the 
(3=2 skeleton simulations in Fig. [9]. Both in the ahistoric context, where ap- 
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proximately half of the nodes change, and with r=19 majority memory, in 
which case only small sets of nodes are changing. 
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Figure 10: Latest snapshots of the changing nodes in one of the j3=2 skeleton 
simulations in Fig. [9]. Color code: Red: 0—s-l, Black: l—>0. Ahistoric patterns 
(left), and patterns with t—19 majority memory (right) . 

Figure [IT] shows the damage spreading in the scenario of Fig. [9]. In this 
scenario, with low connectivity, the depletion in the advance of the damage 
becomes apparent already with r=3 though by T=100 the damage rate also 
reaches the 0.5 level. As expected, r=9 is much more effective in the damage 
control, and, above all, unlimited trailing memory which virtually impedes the 
advance of damage. 
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Figure 11: Damage spreading in the scenario of Fig. [9]. Color code as in Fig. [4]. 



3 Asymptotic levels. 

Figure [12] shows the asymptotic changing rate (left) and damage spreading 
(right) with increasing values of the length of memory r in the three /3-scenarios 
considered before. Asymptotic standing for the mean of then ten last values in 
both cases. The scenario is that of previous sections, i.e., a iV=1000 network 
run up to T = 100 iterations, so that r=100 implies full memory. The changing 
rate notably declines as a result of the initial increase of the length of memory 
(particularly with (5=0.9 ). The lowest changing rate levels are achieved around 
r ~ T/2=50. These levels are not very much altered until a very high length of 
memory is reached. Namely by t > 9T/10=90, the oscillation in the changing 
rate is maintained most of the time-steps, so that the subsequent decay becomes 
increasingly difficult. 



,8=0.9 
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Figure 12: Asymptotic changing rate (left) and damage spreading (right) in 
-ZV=1000 parity /3-skeletons with r-majority memory. 



The effect of increasing the memory length on the damage spreading (Fig. 12 



right), is much more monotone: a seemingly plateau, below the 0.5 level, is 
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reached in the (3=1.0 and /3=2.0 cases, whereas in the highly connected /3=0.9 
context, memory does not qualitatively alter the maximum 0.5 level. 

The effect of memory in bigger networks, does not qualitatively differ from 
that presented here for 7V=1000. The dynamics of both the changing rate and 
damage spreading behave in broad strokes as in the iV=1000 case. Maybe in a 
smoother way in every single simulation, and closer when considering different 
simulations. But smaller networks are more problematic. Thus, Fig 13 shows 



the asymptotic changing rate (left) and damage spreading (right) in the same 
conditions of Fig. [12] except in what refers to the network size : N = 100 instead 
of N = 1000 . The general pattern of the asymptotic changing rate frame in 
Fig.[T3|resembles that of Fig. [12] albeit the higher variability of the changing rate 
in the small network is hidden by the mere presentation of the average of the 
last values. The structure of the damage spreading frame changes dramatically 
compared to that of the bigger network : memory is unable to avoid the full 
damage spread across small networks. Neither in the low connected scenarios, 
not even in that with /?=2.0 . 
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Figure 13: Asymptotic changing rate (left) and damage spreading (right) in 
./V=100 parity /3-skeletons with r-majority memory. 



4 Weighted Memory. 

Majority memory of the last r time-steps demands r extra bits per node to store 
their corresponding state values. To avoid this drawback, past state values can 
be weighted in such a way that only the accumulated memory charge needs to 
be stored. Thus, for example, historic memory can be weighted by applying a 

(T—t) 

geometric discounting process in which the state a\ , obtained t time steps 
before the last round, is actualized to ctV) , a being the memory factor lying 
in the [0,1] interval. This well known mechanism fully takes into account the 
last round (a = 1), and tends to forget the older rounds. 

Every node will be featured by the rounded weighted mean of all its past 
states, so the memory mechanism is implemented in two steps at time-step T : 

(i) The unrounded weighted mean (m) of the states of every node is com- 
puted first : 
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(T), (1) (T)\ 



T-1 

(T) 



t=l 



UJ; 



T-1 

T-t 



n(T) 



Q 

t=l 



(m) Then the trait state s is obtained by comparing the memory charge m 
to the landmark 0.5 (if a £ {0, 1}), assigning the last state in case of an equality 
to this value, so that : 



1 


if 


(T) 

m\ 


> 0.5 




if 


(T) 
m i 

(T) 
m i 


= 0.5 





if 


< 0.5 



The choice of the memory factor a simulates the long-term or remnant 
memory effect: the limit case a = 1 corresponds to a memory with equally 
weighted records (full memory, equivalent to unlimited trailing majority mem- 
ory) , whereas a« 1 intensifies the contribution of the most recent states and 
diminishes the contribution of the more remote states (short-term working mem- 
ory). The choice a = leads to the ahistoric model. This memory implemen- 
tation will be referred to as a-memory. 

This geometric memory mechanism is not holistic but cumulative in its de- 
mand for knowledge of past history: the whole series needs not be known 

(Tl {T} 

to calculate the term ui\ of the memory charge m\ , while to (sequentially) 
calculate w- one can resort to the already calculated u>\ and compute: 
uj\ = abj\ + a I . Consequently, only one number per node (wj) needs 
to be stored. This positive property is accompanied by the drawback of any 
weighted average memory: it computes with real numbers instead of with inte- 
gers as is done with majority memory. 

In the most unbalanced scenario, a\ = ... = a\ T o~\ , it holds thatH: 
m = 0.5 =>• — 2«t + 1 = [i] , where ar holds for the critical value of a below 
which memory has no effect up to time-step T . When T — > oo, the equation [4] 
becomes : —20:00 + 1 = 0, thus, a-memory is not effective if a < 0.5 . 

Figure [14] shows the changing rate with a-memory in parity /3-skeletons of 
one thousand nodes, starting from a single active node. Without memory, the 
changing rate progresses very fast regardless of the /3 value. On the contrary, full 
memory, i.e., a =1.0, dramatically restrains the changing rate. In the /3=2.0 
scenario, thus with low connectivity, even the smallest memory charges, e.g., 
a=0.6, lead the changing rate to extinction. In the more connected networks, 



1 1 a 

m(0,0, ...,0, 1) 



m(l,l,...,l,0) 



2 2 a- 1 

1 a T - 1 1 a T - 1 



2a T + 1 = 



2 a-1 2 a-1 
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low memories such as a=0.6 and a=0.7 only delay the increase of the changing 
rate, that reaches 0.5 in the figure in both the (3—1.0 and (3=0.9 scenarios. With 
a=0.8 memory, the 0.5 level is not reached up to T=100 in the (3—1.0 simulation, 
but it is reached in the (3=0.9 simulation. Finally, with a=0.9 the changing rate 
is notably restrained in both scenarios, though it seems not stabilized as it is 
with full memory. Thus, in the simple context of Fig. [14] full memory exerts a, 
let us say, expected effect, that of maximum dynamic inhibition. The dynamic 
starting from a single active node may be considered as the simplest form of 
damage spreading, measured by the evolution of density. The density curves 



(not shown) mimic those of the changing rate shown in Fig. 14 , including those 
with full memory, agreeing with the expected effect on damage spreading already 
shown in the majority memory scenario. 
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Figure 14: Changing rate with a-memory in parity /3-skeletons from a single 
active node. N=1000 . Red plots correspond to the ahistoric simulations, blue 
plots correspond to full majority memory. 



Figure 15 shows the asymptotic changing rate (left) and damage spreading 
(right) in ./V=1000 parity /3-skeletons with a-memory. The effect of a-memory is 
much more monotone than that of majority memory shown in Figure [12]. Thus, 
as a rule, increasing the memory factor implies restrain in both the changing 
rate and the damage spreading. The changing rate in Fig. 15 sharply peaks 



in the full memory (a=1.0) scenario, according to the the full lengt (r=100) 
majority memory changing rates shown in Fig. [12]. It turns out remarkable that 
damage spreading is restrained only when a reaches high values : > 0.90 in 
the medium connected (3=1.0 scenario, and > 0.85 in the low connected (3=2.0 
scenario. 




Figure 15: Asymptotic changing rate (left) and damage spreading (right) in 
7V=1000 parity /3-skeletons with a-memory. 
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5 Other rules, other contexts 

The results presented up to now deal with the parity rule. This rule is a 
particularly 'active' rule that produces very high changing rates. Neverthe- 
less, most other totalistic rules, i.e., based in the sum of the neighbour states 
5^ = J]] °j T ' ) * tend to led the dynamic to low stable levels of changing rate. 

This feature is particularly visible in majority rule: cr] T+1 -'=l if 2J2i > Ki ; 
cr- T+1 '=0, if 2J2i > Ki where Ki standing for the connectivity of node i . 




Figure 16: Changing rate with the ahistoric majority rule (red), and the ma- 
jority rule with nodes endowed memory of the parity of the last three states 
(blue). (3 =l-skeletons with 10 3 nodes (left), 10 4 nodes (center), and 10 2 nodes 
(right). 

In the simulations presented in Fig. |16| the generic node i preserves its 
last state in case of a tie: a\ —o~\ if = Ki . The dynamics is very 
much sensitive at this respect, so that adopting <7j T+1 — 1 if 2^. = will 

lead the density to 1.0, whereas with cr^ T+1 ^=0 if 2J^. = Ki the density will 
tend to plummet. With the neutral criterion adopted here, there is no reason 



for density drift to any of the extreme values, so that the red lines in Fig. 16 
become stabilized after a short transition period. 

In this scenario, the activity may be supported implementing memory rules 
not of majority type, such as the parity rule acting as memory |18j . Thus for 
example, the parity of the last three states^J as implemented in Fig. 16 . The 
left panel of this figure applies to simulations with N = 10 3 nodes (usually for 
present paper). Simulations with higher number of nodes, e.g., N = 10 4 in 
the central panel, do not qualitatively differ from those with N = 10 3 . They 
present the same general features in their dynamics, maybe with an additional 
characteristic : the evolution of density, changing rate and damage spreading 
is much more similar among them. Opposite, simulations with few nodes, e.g., 
one hundred nodes in the right panel, are much more unforeseeable. Not only in 
what respect to the effect of memory, even the ahistoric dynamics may converge 
to very different densities and changing rates. Thus, for example, in the N = 



10 simulations of the far right panel of Fig. 16 , some ahistoric dynamics are 



led to very low densities of <j=l nodes, e.g., p=0.2, whereas other ones reach 
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high densities, e.g., p—0.8, very distant from the p=0.5 obtained in ahistoric 
simulations with N = 10 3 or N = 10 4 nodes. 

We have studied the effect of memory in several cellular automata and 
Boolean network scenarios [T3J Q31 [H] . The mean connectivity across the eleven 
simulations run here is if =3. 85, K=2.50, and if =6. 64, for /3=1.0, /3=2.0, and 
/3=0.9 respectively. The dynamics studied in this article may be somehow com- 
pared to those simulations in the referenced articles with mean connectivities 
approaching these values. Thus, for example, the way in which a single active 
node propagates activity (a kind of damage spreading) , as shown in the patterns 
at T=100 in Fig. [17]: no restrain in the ahistoric simulation, some dispersion 
with a=0.9, and activity concentrated near the initial active node with full 
memory. 



a=0.0 



a=0.9 



a=1.0 
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Figure 17: Patterns at T=100 with a-memory in the parity , 
ing from a single active node considered in Fig.|14|. 



=l-skeleton start- 



In previous studies we have analyzed the effect of memory in some spa- 
tialized games. Thus, the battle of the sexes [lOj and the prisoner's dilemma 
(PD) [THl E] . It has been concluded in those studies that the consideration of 
previous choices and payoffs in the updating rule, notably boosts cooperation. 
Although the effect of memory on games with players located in /3-skeleton 
nodes is planned for further work, let us point here that the just outlined con- 
clusion, reached in other spatial contexts, may likely also be found when dealing 
with the PD /3-skeletons. Thus for example, starting with a single defector in 
a sea of cooperators, memory will restrain the defection progress, much in the 
way of the restrain of the advance of the active nodes in Fig. 17 . Memory will 



likely have even a stronger effect on the PD dynamics in /3-skeletons with low 
connectivities, i.e., large f3. On the contrary, in highly connected /3-skeletons, 
i.e., small /3, it is likely to be found that memory will be unable to avoid the 
spread of defection propelled via action at distance. 
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6 Potential applications 

What are potential outcomes of our findings?. Proximity graphs, particularly 
relative neighbourhood graphs (/3-skeletons for /3 = 2) are invaluable in simula- 
tion of human-made road networks, see for example the study of Tsukuba central 
district [HJ 22] . The graphs provide a good formal representation of biologi- 
cal transport networks, particularly foraging trails of ants [T] and protoplasmic 
networks of slime mold Physarum polycephalum [3J H] . Our computational ap- 
proaches could be applied toward studies of damage propagation in artificial 
and natural networks, and also used to develop novel techniques for controlling 
space-time dynamics on the networks. 

Majority of the results in the paper deal with parity rule. This rule is im- 
portant because parity rule cellular automata are considered to be low level 
models of self- replication, see e.g. [32 [33], and also, similar rules are used for 
fast simulation of soliton cellular automata, see e.g. [30]. Ultimately the par- 
ity rule can be considered as a good discrete approximation of propagation of 
perturbations in active non-linear discrete media, where a localised traveling 
perturbation (analog of excitation in excitable media) generates other pertur- 
bations. Thus, results of the paper could be, in principle, applied to studying of 
memory-based control of spatio-temporal dynamics of non-linear active media. 

With regards to smile mold, indeed its protoplasmic network is always chang- 
ing, new tubes are formed, some old tubes become abandoned. Even in this case 
- 'a la structure dans le flux' — /3-skeleton is an excellent approximation of 
protoplasmic networks for over 80% of plasmodium foraging span [5] . The parity 
simulated in the paper generates non-trivial patterns of activity, which can be 
interpreted in terms of oscillatorty dynamics of electrical potential or contractil- 
ity of the plasmodium of Physarum polycephalum |39| . By tuning capacity of a 
node memory of /3-skeleton we could shape periodical activity of the graph and 
thus imitate interactions between electrical and oscillatory activity of the slime 
mould, and also incorporate heterogeneous network of discrete bio-chemical os- 
cillators in the discrete models of the plasmodium [7]. 

Yet another, rather more specific domain where our results on /3-skeleton 
with memory can be applied is a computation in assembles of lipid vesicles filled 
with Belousov-Zhabotinsky (BZ) mixture [331 US]- When implementing compu- 
tation in excitable chemical systems [3] we usually encounter an a troublesome 
problem of wave-fragments' instability. A traveling wave-fragment (which rep- 
resents a value of logical variable) rarely preserves its shape for a long time, 
it rather collapses or expands. This problem can be overcome by subdividing 
computing substrate into interconnected compartments — BZ-vesicles — and 
allowing waves to collide one with another only inside the compartments, see 
details in [3 [30] . An ensemble of BZ-vesicles forms a microscopic analog of a 
massively-parallel processor. In computer simulations 9, 20 BZ-vesicles are ar- 
ranged in regular arrays. It is very unlikely that in natural conditions ensembles 
of BZ-vesicles will be regular. When BZ-vesicles, quite possible of different sizes, 
are aggregated into an assemble and tightly packed, they are represented by De- 
launay triangulation. A computation on Delaunay triangulation is implemented 
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in the same manner as on a regular cellular automata [5J. However, when some 
BZ-vesicles start to degrade, and other increase in size average connectivity of 
the vesicle conglomerate decreases. What kind of proximity graph represents 
then the ensemble of BZ-vesicles? 

In 1980 Toussaint demonstrated that MSTCRNGCDT [37]. This hierarchy 
of graphs was later enriched with GG [37] as follows: 

NNG C MST C RNG C GG C DT 

That is when BZ-vesicles bursting the graph representation of the vesicles' 
ensemble will travel downwards in the Toussaint hierarchy: the diagram becomes 
a Gabriel graph GG (/3-skeleton for (3=1) then relative neighbourhood graph 
RNG (/3-skeleton for (3 = 2) then spanning tree MST and then disconnected 
nearest-neighbour graph. Response of each BZ-vesicle to the excitation in the 
neighbouring vesicles depends on the diffusion speed, size of pores between the 
vesicles, concentration of reactants inside the vesicle and overall age of the 
system. Usually, vesicle's reaction in any particular time step is determined 
not only by current dynamics of excitation in the vesicle's neighbourhood but 
also by history of chemical reactions in the vesicle itself. Thus, in computational 
experiments with /3-skeletons with memory we discover, at least at the level of 
discrete yet fine-grained models, how activity can propagate in a ensemble of 
'ageing' BZ-vesicles and what methods of the activity control can be employed. 

7 Conclusion 

As a rule, endowing nodes in beta-skeleton parity automata with majority mem- 
ory induces a moderation in the rate of changing nodes and in the damage 
spreading, albeit in the latter case memory turns out ineffective in the control 
of the damage with low memory charge and in networks with high connectivity, 
i.e., with low (3 . In a simulation up to T time-steps, the maximum depleting 
effect of memory on the changing rate is achived with a memory length r ~T/2. 
With very high memory length, the oscillatory-like effect induced by memory 
in the changing rate is maintained too long to allow a effective depletion in the 
changing rate. Thats is so even in low-connected networks, i.e., with large (3 . 
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